J Arid Land 
https://doi.org/10.1007/s40333-018-0017-y 


Canopy cover or remotely sensed vegetation index, 
explanatory variables of above-ground biomass in 
an arid rangeland, Iran 


Fatemeh PORDEL!, Ataollah EBRAHIMI", Zahra AZIZI? 


! Department of Range and Watershed Management, Shahrekord University, Shahrekord 8818634141, Iran; 
2 Department of RS-GIS, Science and Research Branch, Islamic Azad University, Tehra 1477893855, Iran 


Abstract: Estimation of above-ground biomass is vital for understanding ecological processes. Since 
direct measurement of above-ground biomass is destructive, time consuming and labor intensive, canopy 
cover can be considered as a predictor if a significant correlation between the two variables exists. In this 
study, relationship between canopy cover and above-ground biomass was investigated by a general linear 
regression model. To do so, canopy cover and above-ground biomass were measured at 5 sub-life forms 
(defined as life forms grouped in the same height classes) using 380 quadrats, which is 
systematic-randomly laid out along a 10-km transect, during four sampling periods (May, June, August, and 
September) in an arid rangeland of Marjan, Iran. To reveal whether obtained canopy cover and 
above-ground biomass of different sampling periods can be lumped together or not, we applied a general 
linear model (GLM). In this model, above-ground biomass was considered as a dependent or response 
vatiable, canopy cover as an independent covariate or predictor factor and sub-life forms as well as 
sampling periods as fixed factors. Moreover, we compared the estimated above-ground biomass derived 
from remotely sensed images of Landsat-8 using NDVI (normalized difference vegetation index), after 
finding the best regression line between predictor (measured canopy cover in the field) and response 
vatiable (above-ground biomass) to test the robustness of the induced model. Results show that 
above-ground biomass (response variable) of all vegetative forms and periods can be accurately predicted 
by canopy cover (predictor), although sub-life forms and sampling periods significantly affect the results. 
The best regression fit was found for short forbs in September and shrubs in May, June and August with 
R? values of 0.96, 0.93 and 0.91, respectively, whilst the least significant was found for short grasses in 
June, tall grasses in August and tall forbs in June with R? values of 0.71, 0.73 and 0.75, respectively. Even 
though the estimated above-ground biomass by NDVI is also convincing (R?=0.57), the canopy cover is a 
more reliable predictor of above-ground biomass due to the higher R? values (from 0.75 to 0.96). We 
conclude that canopy cover can be regarded as a reliable predictor of above-ground biomass if sub-life 
forms and sampling periods (during growing season) are taken into account. Since, (1) plant canopy cover 
is not distinguishable by remotely sensed images at the sub-life form level, especially in sparse vegetation 
of arid and semi-arid regions, and (2) remotely sensed-based prediction of above-ground biomass shows a 
less significant relationship (R?=0.57) than that of canopy cover (R? ranging from 0.75 to 0.96), which 
suggests estimating of plant biomass by canopy cover instead of cut and weighting method is highly 
recommended. Furthermore, this fast, nondestructive and robust method that does not endanger rare 
species, gives a trustworthy prediction of above-ground biomass in arid rangelands. 
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1 Introduction 


Above-ground biomass is the amount of vegetative parts produced or energy stored per unit of 
area (Sala and Austin, 2000; Flombaum and Sala, 2007). The study of above-ground biomass of 
rangeland has been a central theme in ecological studies for evaluating their spatial and temporal 
variations controlled by biotic and abiotic factors (Olofsson et al., 2001; Radloff and Mucina, 
2007), and addressed by a range of questions, from livestock forage availability (Di Bella et al., 
2005) to global carbon cycle balance (Sala and Austin, 2000), as well as the prediction of global 
climate change (Li et al., 2014). Therefore, accurate estimation of above-ground biomass is 
critical for the management of rangeland, especially when the effects of management decisions 
need to be predicted (Ebrahimi et al., 2008) or dealt with sustainable management (Snyman, 
1998). Different methods have been developed for measurement of above-ground biomass but 
most common techniques for biomass estimation consist in clipping and weighing a number of 
representative plots or individuals (Montés, 2009), which is time-consuming and costly, even 
non-reproducibility over time (Bonham, 2013), and the species being on the brink of extinction 
are all among the critics against this technique from the economic and scientific point of view 
(Montés, 2009). 

Plant destruction is an important risk, particularly in some ecosystems where the mean density 
of several preferred species is low (Guevara et al., 2002). This issue is more concerned in arid and 
semi-arid rangelands where frequently monitoring of vegetation canopy cover via clip and 
weighting method might cause extinction of rare species in one hand and plant canopy cover 
doesn't show a strong correlation by derived vegetation indices from remotely sensed images due 
to sparsity of canopy cover occurrence on the other hand. Hence, we need an accurate, rapid, and 
reliable nondestructive method of estimating above-ground biomass. 

Canopy cover is defined by the Range Inventory Standardization Committee (RISC) as the 
percentage of ground covered by a vertical projection of the outermost perimeter of the natural 
spread of foliage of plants, while foliage cover projected aerial parts of vegetation onto the 
ground (Anderson, 1986), and only takes into account intercepts of leaves or branches. 

Non-destructive methods have been developed specifically for shrubberies (Etienne, 1989; 
Montes et al., 2004; Frank et al., 2005) and grasslands (Thursby et al., 2002; Tackenberg, 2007). 
For example, canopy cover as an indirect parameter was used to estimate Atriplex verrciferum M. 
B. and Salsola dendroides Pall. biomass at Gherekhlar region in Marand, Iran (Mokhtarias! and 
Mesdaghi, 2008). Yang et al. (2017) also successfully developed some allometric models to 
predict biomass of 12 shrub species in Chinese desert grassland. The results of investigating the 
relationship between above-ground biomass and indirect variables such as canopy cover, foliage 
cover, and basal cover in three habitats of grass-shrubland, grassland, and shrubland revealed that 
both canopy cover and foliage cover have a significant relationship with production in most 
species (Arzani et al., 2008), if morphological characteristics of the species are taken into account 
(Adler et al., 2004; N'avar et al., 2004; Foroughbakhch et al., 2005). Although considerable 
research has been devoted to the relationship between plant production and allometric equation of 
vegetation parts such as canopy cover ( Foroughbakhch et al., 2005; Arzani et al., 2008; Arzani et 
al., 2011; Bonham, 2013), rather less attention has been paid to the effects of life-forms and 
sampling dates on this relationship. Many researchers believe equations set based on the 
relationship between above-ground biomass and canopy cover gained at a given time could not be 
generalized to estimate the above-ground biomass in other times (Hughes et al., 1987; Arzani et 
al., 2011). Hence, these changes in canopy cover might influence its relation with above-ground 
biomass. Therefore, the present paper aims at addressing partially this complex problem in an arid 
rangeland. 

In addition to ground measurement of vegetation characteristics, remotely sensed data provide 
opportunity to monitor large surfaces, regularly with comparatively reasonable times and costs 
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(Lu et al., 2004; Caprioli et al., 2006). Multi-temporal remotely sensed images can be used to 
estimate above-ground biomass and its fluctuations from vegetation indices (VIs). Amongst 
others, Landsat-8 satellite, provides valuable images for investigating land surface features (Ke et 
al., 2015), more specifically, above-ground biomass. Most researchers calculate NDVI 
(normalized difference vegetation index) probably the most used VIs in rangeland and forest 
management (Brinkmann et al., 2011; Migliavacca et al., 2011; Zhu and Liu, 2015). Several 
researchers have shown that the relationship between vegetation and environmental parameters 
with VIs can be expressed as an empirical model (Du Plessis, 1999; Drake et al., 2003; 
Brinkmann et al., 2011; Ji et al., 2012; Porter et al., 2014; Zhu and Liu, 2015). Glenn et al. (2016) 
predicted accurately above-ground biomass with remotely sensed data in western USA at some 
life-form and total above-ground biomass levels. Pordel (2015) established some models to depict 
spatio-temporal changes of above-ground biomass using VIs derived and extrapolated from 
Landsat-8 images in Mrajan rangeland, Iran. Comparing the measurement of above-ground 
biomass in the field with the satellite-based derived above-ground biomass will also reveals pros 
and cons of these valuable methods that will be also partially deals with in this research. In 
summary, a central issue in the study of relationship between above-ground biomass and canopy 
cover or VIs is studying the effects of life-forms and sampling periods that may influence this 
relationship due to changes in plant composition or phenological stages, an important issue that 
received few attentions so far. Therefore, we can summarize the objectives of the preset study as: 
(1) examining whether the relationship between canopy cover (predictor) and above-ground 
biomass (response) is influenced by sampling periods and sub-life forms or not; (2) fitting the 
best prediction equation between above-ground biomass using indirect indicator of canopy cover 
at different sub-life forms and sampling periods; and (3) predicting and comparing above-ground 
biomass derived from NDVI with the evaluated above-ground biomass by means of indirect 
explanatory variable of canopy cover. 


2 Materials and methods 


2.1 Study area 


This study was conducted in the Marjan rangelands of Boroujen County, Iran, which is located at 
31°58'51"N to 32°08'29'N and 51°16'50"E to 51°24'38”E, with a mean elevation of 2200 m a.s.l. 
(Fig. la). According to the Koppens' classification method (Wang and Overland, 2004), this area 
has a temperate and cold climate along with warm and dry summers. Long-term mean annual 
precipitation (1988-2013) was 250 mm, falling mainly in fall and winter with the highest monthly 
mean precipitation occurring in December (48.8 mm) and January (48.1 mm) and the lowest in 
Augusts (0.3 mm) and September (0.6 mm; Fig. 1b). The study area is approximately 7830 hm?. 
Some dominant plant species, along with their family and sub-life forms are presented in Table 1. 


2.2 Methods 


2.2.1 Field measurements of canopy cover and above-ground biomass 


Canopy cover (i.e., outermost perimeter of the plants) and above-ground biomass were estimated 
using 95 quadrats at 19 sampling sites (each contains 5 quadrates, with 1 centroid quadrat and 4 
quadrates positioned on the four directional corners with 5-m intervals from the central quadrat), 
the so called sampling nodes, that systematic-randomly distributed along a 10-km transect (the 
first node was selected systematically but the rests were randomly distributed along the transect) 
(Fig. 2). The aforementioned measurements were repeated for 4 sampling periods during spring 
and summer seasons in 2015 (i.e., May, June, August, and September), which made 380 sampling 
quadrates in total. Following the method of Tahmasebi et al. (2012), a 2 mx2 m quadrate size was 
chosen for estimating canopy cover and above-ground biomass. Canopy cover and above-ground 
biomass were measured at 5 sub-life form levels, defined as the morphologically similar species 
from the same vegetation stratum (height class) (Ebrahimi et al., 2008). 
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Fig.1 Sampling sites (a) and monthly mean temperature and precipitation (b) of the study area 


Table 1 Dominant plant species, their genera, families, and sub-life forms in the Marjan rangeland 


Family Genus Species Sub-life form 
Poaceae Bromus Bromus tomentellus Boiss. Tall grass 
Poaceae Festuca Festuca ovina L. Tall grass 
Poaceae Stipa Stipa hohenackeriana Trin. & Rupr. Tall grass 
Poaceae Taeniatherum Taeniatherum crinitum Schreb. Short grass 
Poaceae Bromus Bromus tectorum L. Short grass 
Poaceae Bromus Bromus danthoniae Trin. Short grass 

Chenopodiacea Noaea Noaea mucronata (Forssk.) Shrub 
Astraceae Scariola Scariola orientalis (Boiss.) Sojak Shrub 
Papilionaceae Astragalus Astragalus verus Olivier Shrub 
Papilionaceae Astragalus Astragalus cephalanthus DC. Shrub 
Papilionaceae Astragalus Astragalus pinetorum Boiss. Short forb 
Brassicaceae Alyssum Alyssum marginatum Steud. ex Boiss. Short forb 
Papilionace Onobrychis Onobrichis gaubae Bornm. Short forb 
Lamiaceae Phlomis Phlomis persica Boiss. Tall forb 
Lamiaceae Phlomis Phlomis olivieri Benth. Tall forb 
Lamiaceae Stachys Stachys pilifera Benth. Tall forb 


Cut and weight method was used to measure above-ground biomass (Bonham, 2013) at the 
sub-life form level. Above-ground biomass was cut in quadrats, green and photosynthetic leaves 
and twigs were separated from dried parts at each vegetative period. Green above-ground biomass 
was oven-dried at 65°C and finally weighed. A gridded quadrat frame of 2 mx2 m was used to 
accurately measure canopy cover at the sub-life form level (Shiyomi and Yoshimura, 2000). Each 
sides of quadrate were divided into 100 equal parts of 20 cmx20 cm as a mesh or gridded 
quadrate to be more precise in estimating of canopy cover (Bonham, 2013). For each sampling 
period only green canopy cover was measured. 

2.2.2 Data analysis of modelling the relationship between canopy cover and above-ground 
biomass 


The process of data analysis is summarized as underneath stages. At the first stage, normal 
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Fig.2 Sampling design for testing the relationship between canopy cover and above-ground biomass. Schematic 
representation shows 3 out of the 19 sampling sites (with 400-1000 m intervals) along a 10-km transect. Each 
sampling site contains 5 quadrats of 2 mx2 m dimensions arranged with 1 centroid quadrat and the other 4 
quadrats positioned 5-m away at four directional corners. 


distribution of response variable (above-ground biomass) was investigated. We applied a general 
linear model (GLM) to test if a significant relationship between above-ground biomass and 
canopy cover exists and whether this relation is affected by life form and sampling date or not. In 
the GLM, measured above-ground biomass was regarded as the dependent variable (response 
variable), canopy cover as the independent covariate (predictor or explanatory variable) and 
sub-life form and sampling periods as fixed factors. 

After testing the relationship between response and predictor variables, we fit the best 
regression model between dependent or response (above-ground biomass) and independent or 
predictor variable (canopy cover) at different sub-life forms and sampling dates. We developed 
linear predictive models for dependent variable of above-ground biomass based on independent 
variable of canopy cover using a set of 60% of randomly selected data out of all the ground 
sampled data (destructively sampled quadrats). 

After processing the regression models, we investigated the following items that are the 
primary default in regression models to distinguish robustness of each predictive model: (1) 
Normal Distribution of Residuals Frequencies (NDRF) by normal P-P plot of standardized 
residuals; (2) testing a non-correlation between the residuals and the Durbin-Watson test (Durbin 
and Watson, 1951); (3) depicting the Stability in Variances of the Residuals (SVR) in a scatter plot 
where standardized residuals showed in Y axis and the standard predicted (estimated response) 
values in X axis to detect non-linearity and inequality error of variances; and (4) looking over the 
robustness of Linear Regression (LR) relationship between Y (dependent variable) and X 
(independent variable) by R? and r. 

2.2.3 Modelling the relationship between NDVI and above-ground biomass 


After selecting the most cloud-free and nearest acquiesced images of Landsat-8 satellite to the 
sampling dates over the study area (Table 2), we applied FLAASH atmospheric correction on the 
images and calculated NDVI from the resulted images based on Equation 1. 
_ NIR -RED (1) 

NIR + RED’ 
where NIR is the reflectance in the near-infrared band; and RED is the reflectance in red band 
(Tucker, 1979). 

The values of NDVI at sampling nodes were extracted and regressed against corresponding on 
field measured above-ground biomass of May to September (all sampling dates), where 
above-ground biomass was considered as a response variable and NDVI values as a predictor or 
independent variable. The primary defaults of regression models (mentioned in Section 2.2.2) 
were investigated to distinguish robustness of the predictive models. 

Finally, a pair-sampled t-test was performed between derived above-ground biomass from 
NDVI values and on field measured above-ground biomass to determine if a statistically 
significant difference exists between actual and estimated above-ground biomass values (in 40% 
of the data that regarded as the test) for both methods of remotely sensed- and canopy 


NDVI 
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cover-based methods. All of statistical analysis was performed with SPSS software (version 
22.00). 


Table 2 Landsat 8 images used for calculating vegetation indices and investigating its relation with 
above-ground biomass 


Designation Date (yyyy-mm-dd) Cloud cover (%) 
A 2014-05-24 8 
B 2014-06-25 3 
C 2014-08-02 0 
D 2014-09-19 0 
3 Results 


3.1 Field measurement 


Canopy cover and above-ground biomass in all sub-life forms and sampling periods are presented 
in Figure 3. In total, averages of above-ground biomass and canopy cover decreased from 
sampling period of May to September. Shrub had the highest averages of canopy cover and 
above-ground biomass in all sampling periods, but canopy cover of the short forbs was only 
larger than that of the shrubs in May (Fig. 3). Average canopy covers of all sub-life forms were 
noticeably decreased from May to September unless the shrubs that slowly increased from May to 
June then slightly decreased from June to August and finally decreased by 1.5% from August to 
September. In this respect, short forbs and short grasses showed a more sever decreasing from 
May to September in comparison to the other sub-life forms. When above-ground biomass was 
concerned, the same pattern was observed, however, increasing of shrub biomass from May to 
June was more outstanding than that of canopy cover. Moreover, the biomass of shrubs was 
highest throughout all sampling periods whilst its canopy cover was slightly lower in May than 
that of tall grasses. Generally, we can conclude that shrubs are dominant in terms of both canopy 
cover and above-ground biomass in the study area. Furthermore, short grasses (i.e., Taeniatherium 
crinitum Schreb, Bromus tectorum L., and Bromus danthoniae Trin.) were absent in August and 
September, whilst tall grasses was only occurred in September (Fig. 3). Overall, above-ground 
biomass and canopy cover of shrubs were considerably higher than that of other sub-life forms 
during sampling periods. Tall forbs and short grasses showed the lowest values in this period (Fig. 
4). 
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Fig. 3 Averages of above-ground biomass (a) and canopy cover (b) from field measurement in all sub-life forms 
and sampling periods 


3.2 GLM outputs 


Results of GLM analysis of above-ground biomass and canopy cover as well as the main effects 
of sub-life forms and sampling periods and their interactions are shown in Table 3. We found a 
significant correlation (R?=0.81) between above-ground biomass (response variable) and canopy 
cover (predictor variables) for all sub-life forms and sampling periods. Sampling periods and 
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Fig. 4 Averages of total above-ground biomass (a) and canopy cover (b) during sampling periods at different 
sub-lifeform levels 


sub-life forms as well as interaction of these two factors showed a significant effect on estimation 
of above-ground biomass (P<0.05). A significant difference (P<0.05) also found in the intercept 
of the linear regression of the models (Table 3). 


Table 3 ANOVA of the general linear model (GLM) applied on dependent (above-ground biomass) and 
independent covariate variable (canopy cover) with main effects of sub-life form and sampling period (fixed 
factors) and their interactions 


Source Sum of square df Mean square F P 
Corrected model 1,125,513.40 17 66,206.60 171.80 0.001 
Intercept 3879.01 1 3879.01 10.00 0.002 
Canopy cover 735,494.30 1 735,494.30 1905.40 0.001 
Life form 43,303.90 4 10,825.90 28.04 0.001 
Sampling period 14,121.30 3 4707.10 12.10 0.001 
Life formxSampling period 17,085.70 9 1898.40 4.90 0.001 
Error 255,910.30 663 385.90 
Total 2,035,461.60 681 
Corrected total 1,381,423.70 680 


3.3 Modelling the relationship between canopy cover and above-ground biomass 


In general, the results of linear regression in all sub-life forms and sampling periods indicated a 
high R? value (i.e., 0.71 for short grasses and 0.96 for short forbs; Table 4). Overall, shrubs had 
the most significant relation (mean of R7=0.91), short grasses and tall grasses showed the least 
significant relation (mean of R?=0.78), and tall forbs and short forbs lay in between with the mean 
R? values of 0.80 and 0.84, respectively. R? changes from sampling periods of May to September. 
In total, September, May, August, and June had the highest to the lowest values of mean R? from 
0.94, 0.85, 0.82, to 0.80, respectively. As indicated in Table 3, all of the primary defaults of 
regression models were met. 


3.4 Modelling the relationship between NDVI and above-ground biomass 


Table 5 shows a linear regression between above-ground biomass (response variable) and NDVI 
(predictor variable). A significant relationship is also found between above-ground biomass and 
NDVI (P<0.05) with a R? value of 0.57. As indicated in Table 5, all of the primary default of 
regression model in this method was also met. 


3.5 Results of validity test of the models 


The results of estimated above-ground biomass from developed prediction models versus field 
measurements of above-ground biomass are shown in Figure 5. Results reveal that no significant 
(P>0.05) difference is found between predicted above-ground biomass derived from developed 
models of canopy cover and field measurements of above-ground biomass for all sub-life forms 
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and sampling periods in testing samples (Figs. 5a—e). Similarly, no significant (P>0.05) difference 
is found between predicted total above-ground biomass derived from NDVI and measured 
above-ground biomass in the field (Fig. 5f). 


Table 4 Developed predictive regression models between canopy cover (independent) and above-ground 
biomass (dependent) variables as well as the results of primary defaults in regression models 


Sampling date Sub-life form Equation R NDRF VRS Durbin-Watson LR of 
Tall grasses Y=6.3+8.03X 0.84 v v 1.50 v 
Tall forbs Y=4.33+4.39X 0.77 v v 2.20 v 

May Short forbs Y= —0.29+8.81X 0.88 v v 1.40 v 0.85 
Shrubs Y= —0.29+11.56X 0.93 v v 2.00 v 
Short grasses Y=4.07+3.97X 0.83 v v 2.20 v 
Tall grasses Y=4.5+7.52X 0.78 v v 1.60 v 
Tall forbs Y= —0.78+9.01X 0.75 v v 1.65 v 

June Short forbs Y=6.25+3.97X 0.77 v v 2.55 v 0.78 
Shrubs Y= —4.6+18.42X 0.91 v v 2.00 v 
Short grasses Y=1.43+3.14x 0.71 v v 1.70 v 
Tall grasses Y= —0.6+12.14X 0.73 v v 1.60 v 

mad Tall forbs Y= 1.774+4.17X 0.80 v v 2.40 v öz 
Short forbs Y=1.51+3.53X 0.75 v v 2.00 v 
Shrubs Y= —1.86+17.92X 0.91 v v 1.70 v 

Tall forbs Y=-1.7+12.81X 0.88 v v 1.82 v 0.90 

September Short forbs Y= 1.39+3.26X 0.96 v v 1.70 v 
Shrubs Y= —0.27+16.72X 0.87 v v 2.23 v 


Note: Y, above-ground biomass X, canopy cover; NDRF, Normal Distribution of Residuals Frequencies; SRV, Stability of Residuals 
Variances; Non-correlation between the residuals: 1.50<Durbin-Watson<2.50; LR, linear relationship between dependent and 
independent variables. v indicates that primary default of regression model was met. It should be noted that some of the sub-life 
forms were absent in some sampling periods. 


Table 5 Linear model between NDVI and total above-ground biomass 
Equation R? NDRF SVR Durbin-Watson LR Sig 


Y= —47.64+1273.06xNDVI 0.57 v v v v 0.00 


Note: NDRF, Normal Distribution of Residuals Frequencies; SVR, Stability of Residuals Variances; Non-correlation between the 
residuals: 1.50<Durbin-Watson<2.50; LR, linear relationship dependent and independent variables. v indicates that primary default 
of regression model was met. 


3.6 Above-ground biomass images derived from NDVI and measured values 


The above-ground biomass maps obtained from the best fit of linear models between NDVI and 
total above-ground biomass at different sampling periods are shown in Figure 6. As indicated in 
these images, NDVI values as well as above-ground biomass generally declined from May to 
September. The highest values of above-ground biomass observed in May ranged from 300 to 600 
kg/hm? (Fig. 6a), while the lowest values observed in September were between 20 and 140 
kg/hm? (Fig. 6e). 


4 Discussion 


In this study, we investigated that whether canopy cover could be an appropriate explanatory 
variable for estimating above-ground biomass in an arid steppe rangeland. The results were 
compared with the remotely sensed data. As shown in Figure 4, the results from the measurements 
in the field indicated that the averages of canopy cover and the above-ground biomass decreased 
from 15.70% to 2.85% and from 330.3 to 97.5 kg/hm? from May to September, respectively. 
Declining vegetation from the end of May in the Mediterranean arid steppe regions is a normal 
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Fig. 5 Predicted above-ground biomass and total above-ground biomass from developed models versus 
measured above-ground biomass in the field at different sub-life forms and sampling periods using testing data. 
There is no significant difference between predicted and measured values (P>0.05). 


phenomenon due to the lack of rainfall during summer season (Fig. 1b). In the study area, canopy 
cover and above-ground biomass of shrubs are dominated in comparison to the other sub-life 
forms (Figs. 3 and 4). The fluctuation of this life form in above-ground biomass and more 
specifically in canopy cover is lower than those of any other sub-life forms. Namely, other 
sub-life forms significantly decreased in both canopy cover and above-ground biomass from May 
to September (Fig. 3). In contradict to that, the shrub life form tended to increase from May to 
June and August and finally decreased to September. The growth season of shrub species (e.g., 
Astragalus verus Olivier and Astragalus cephalanthus DC.) starts with a significant delay 
compared to other life forms. The growth season of short forbs (mostly ephemerals) grows up 
earlier than those of tall grasses. Intrinsically, short forbs (e.g., Astragalus pinetorum Boiss., 
Astragalus effuses Bunge, Allysum marginatum Steud. ex Boiss., and Onobrychis gaubae Bornm 
in Table 1) are mostly replaced by other sub-life forms during the growth season, especially 
shrubs and tall grasses that had more stability in above-ground biomass. Although, short forbs 
dominated the first sampling course, still they seriously declined in above-ground biomass and 
canopy cover in comparison to the shrubs and tall grasses during the second and third sampling 
periods. Moreover, we found a few short forbs in the last sampling periods. However, tall forbs 
and short grasses were absent in the last sampling period (Figs. 3 and 4). 

As shown in Table 3, a significant relation (P<0.05) was found between above-ground biomass 
(response variable) and canopy cover (predictor variable), suggesting that canopy cover is an 
instructive explanatory variable for estimating above-ground biomass. Regardless of 
aforementioned fluctuation in biomass and canopy cover of different sub-life forms which is a 
commonplace phenomenon in arid and semi-arid rangeland of the Mediterranean type, the result 
of Table 3 indicates that above-ground biomass can accurately be predicted by canopy cover. 
Montes et al. (2004), Flombaum and Sala (2007) and Arzani et al. (2008) also found that canopy 
cover is an appropriate proxy for biomass estimation. Results from Table 4 and Figure 5 also 
demonstrate the hypothesize that canopy cover properly describes the above-ground biomass 
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Fig. 6 Above-ground biomass (AGB) maps at different sampling periods derived from calculated NDVI. The 
white area in the study area is the other land uses. 


in arid rangeland where canopy cover is distinct and sparse. Thereby, this indirect method of 
estimating plant biomass that is based on an easier-to-measure attribute of canopy cover is vital 
for regularly tuning of grazing pressure. This method is quite faster and less labor intensive to 
estimate than destructive method of clipping and weighting that endangers rare species. 

Although a high correlation was found between above-ground biomass and canopy cover, the 
data of different sampling dates and sub-life forms should not be lumped together to make one 
predictive model due to a significant effect of sampling dates and sub-life forms, as well as, their 
interaction on estimation of above-ground biomass (P<0.05) by canopy cover (Table 3). This is 
quite dangerous in different growth stages, wherein the height and volume changes are 
perceptible. Therefore, we emphasize the effects of both sampling periods and vegetative forms 
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on the relationship between above-ground biomass and canopy cover. Some literature also 
highlighted the effects of climate change and grazing impacts on the above-ground biomass and 
canopy cover relationships (Hughes et al., 1987; N'avar et al., 2004, Foroughbakhch et al., 2005; 
Arzani et al., 2008; Arzani et al., 2011). We do not recommend the use of predictive models of 
mid-growth stages (June) for other stages (see Tables 2 and 3). Moreover, a significant difference 
(P<0.05) was also found in the intercept of the GLM (Table 2), enforcing us to use different 
prediction models for sampling dates and sub-life forms. The significant difference in intercept of 
GLM might originate from intrinsic variation in above-ground biomass and canopy cover of 
different growth stages (sampling dates) of various sub-life forms. 

Even though the correlation between canopy cover and above-ground biomass is significant in 
all sampling periods (R? ranging from 0.71 to 0.96), still shrubs tended to show a more stable and 
higher relationship between its canopy cover and above-ground biomass (mean of R7=0.90). This 
is particularly when they are being compared to short grasses that showed the least significant 
relationship (mean of R?=0.77) (Table 4). This might be due to a more distinguished canopy cover, 
and a higher volume of shrubby species and a more stable and homogenous outer-form of shrub 
sub-life form in contrast to the others. According to our findings, when the plants have spent their 
maximum growth and greenness, progressing towards dehydration; the correlation between 
canopy cover and above-ground biomass slightly reduced (Table 4). The increasing trend of 
significance that is found from sampling period of August to September might be due to 
increasing shrub life form from early vegetative growth to maturing as well as diminishing short 
grasses that are hard to be predicted. The high correlation at the end of growth season in some 
sub-life forms, such as tall forbs, might be due to the presence and regrowth of new species that 
are green and pass their early and middle periods of growth stages. Other research also found that 
the species-specific allometric models fit the data well to predict total above-ground biomass for 
each species, but the above-ground biomass of tree shaped species was significantly predicted by 
a single predictor. Moreover, shrubs with a tall stem and an umbrella-like canopy structure (e.g., 
Acacia mellifera (Vahl) Benth) were most accurately predicted by use of a combination of both 
circumference of the stem at ankle height and canopy volume (Hasen Yusuf et al., 2013). A higher 
correlation found for shrubby species in our finding is in line with Hasen Yusuf et al. (2013), who 
approving a more reliable predicting of shrub species than the others. 

In this study, we found no significant difference between the estimation of above-ground 
biomass by using NDVI data and ground sampling in the field (Fig. 6). Meanwhile, a relatively 
poor relation between above-ground biomass and NDVI at the sub-life form level was found. 
Ghorbani et al. (2017) also found a significant but no strong relationship (R7=0.23, 0.29 and 0.153 
for grasses, forb and shrubs, respectively) between above-ground biomass and NDVI in Sabalan 
rangelands of Iran. This might be due to the effect of factors such as soil reflectance in arid and 
semi-arid regions on the received signals by the sensors of satellites. 

Since no significant difference between above-ground biomass estimates using canopy cover 
and field measurements at different sub-life forms and sampling periods (Fig. 5) was found, it can 
be concluded that canopy cover is an appropriate explanatory variable for predicting 
above-ground biomass. Still the sampling periods and sub-life forms should be taken into account 
due to morphological variations during different phenological stages and instinct variation of 
sub-life forms canopy cover. Since the term canopy is referring to the extent of the outer layer of 
leaves of an individual plant species or a group of plants that normally blocks light to the ground, 
its structure, organization, or spatial arrangement (three-dimensions) might differ on phenological 
stages from early vegetative growth to maturity. This might be the main reason of significant 
effect of sub-life form and sampling date upon the relationship between predictor and response 
variables (Table 3). Recent research suggests that canopy cover is a reliable predictor variable and 
an appropriate alternative indicator for estimating above-ground biomass at the sub-life form level 
(Ebrahimi, 2017); this potential generally can help to estimate rangelands biomass and grazing 
capacity as a non-destructive method (Gholami baghi et al., 2013). 

The results revealed that despite the capability of remotely sensed images for estimating 
above-ground biomass by means of NDVI (Figs. 5 and 6; Tables 4 and 5), canopy cover is a more 


JOURNAL OF ARID LAND 


straightforward and reliable indicator for predication of above-ground biomass not only at the 
community level but also at the sub-life form level (R? ranging from 0.75 to 0.96). Therefore, we 
recommend the use of this nondestructive, fast, and reliable method for estimating above-ground 
biomass in arid and semi-arid rangelands. 


5 Conclusions 


Due to spare distribution of vegetation in arid and semi-arid rangelands, estimation of 
above-ground biomass by remotely sensed data is not feasible, however, using canopy cover to 
estimate above-ground biomass especially the shrubby species is highly recommended for its 
non-destructive. According to our research, there is a strong relationship between canopy cover 
and above-ground biomass. Therefore, estimation of above-ground biomass is accurately possible 
using canopy cover at the sub-life form level. However, given the effects of growth season and 
sampling periods on this relationship, developing different models to predict above-ground 
biomass based on explanatory variable of canopy cover is needed. 
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